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Abstract 

We study scattering properties of a PT-symmetric square well poten- 
tial with real depth larger than the threshold of particle-antiparticle pair 
production as the time component of a vector potential in the Dirac equa- 
tion. Spontaneous pair production inside the well becomes tiny beyond 
the strength at which discrete bound states with real energies disappear, 
consistently with a spontaneous breakdown of PT symmetry. 

1 Introduction 

The study of PT-symmetric potentials originated in the seminal papers [I], [5] 
by Bender and coworkers, dedicated to the analysis of spectra of non-Hermitian 
Hamiltonians of anharmonic oscillators, which turn out to be entirely real on 
condition that exact PT symmetry holds, i.e., the Hamiltonian, H, commutes 
with the PT operator and all eigenfunctions of H are also eigenstates of VT. In 
case of exact VT symmetry, it is possible to formulate an equivalent quantum 
mechanical description by defining a new metric operator in the representa- 
tion space, a concept already established [3] before the study of the properties 
of 'PT-symmetric potentials and generalized later to the definition of pseudo- 
Hermitian quantum mechanics [4]. If the condition that the eigenfunctions of the 
VT invariant Hamiltonian are also eigenfunctions of VT is relaxed, the VT sym- 
metry is spontaneously broken and complex eigenvalues appear in the spectrum 
of H. Since then many authors have examined bound state problems with 
PT-symmetric potentials in non relativistic quantum mechanics, while rela- 
tively few studies are dedicated to relativistic models, or even to PT-symmetric 
quantum field theories [SI [B]. The present work intends to be a contribution 
to the study of the VT invariant Dirac equation, with emphasis on scattering 
aspects, within the framework of standard relativistic quantum mechanics. PT- 
symmetric potentials will be treated phenomenologically as effective potentials, 
without attempting to formulate an equivalent Hcrmitian theory [6]. 
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Solvable models in the Dirac equation are most easily constructed in (1+1) 
space-time dimensions and the related 'PT-symmetric potentials may have dif- 
ferent behaviour under Lorentz transformations, and are classified as vectors, 
pseudovectors, scalars, or pseudoscalars. In the cases of scalar and pseudoscalar 
potentials, affecting particles and antiparticles in the same way, the relevant hid- 
den symmetry of the Dirac equation can be classified as pseudosupersymmetry[7]. 
The analysis of vector potentials has been focused mainly on bound states, like 
in the generalized Hulten potential of Ref . [5] , or in the logarithmic derivative of 
a suitable position-dependent effective mass of Ref. [3] , while scattering aspects 
have not been examined in detail. Among the latter, an interesting peculiarity 
of strong potentials, with \V\ > 2m, the threshold for production of a particle- 
antiparticle pair, with m the particle mass, is the possibility that bound states 
merge with the negative-energy continuum of scattering states, thus appearing 
as transmission resonances at negative energies inside the potential well, which 
becomes overcritical with respect to spontaneous creation of particle- antiparticle 
pairs. 

Such overcritical vector potentials are well studied in the Hermitian case in 
(1+1) dimensions, from the simplest example of the square wellflO]. to the cusp 
potential [TT], [T2], or the Woods-Saxon potential[T3]. In the present work, we 
study scattering aspects of overcritical VT -symmetric potentials, in particular 
how overcriticality is interrelated with spontaneous breakdown of VT symme- 
try. It is worthwhile to remark from the very beginning that our framework 
is standard quantum mechanics with complex potentials, phenomenologically 
treated as effective potentials. For the sake of simplicity, the analysis is carried 
out in detail for the PT-symmetric square well, taken as the time component 
of a vector potential. 

In this short note, we shall not enter into discussion of construction of pos- 
itive definite norms via a linear operator, C, commuting with the VT operator 
and the Hamiltonian, H, introduced by Bender and collaborators in nonrel- 
ativistic quantum mechanics [Tl] . [T5] . [Tfi] and later extended to quantum field 
theory, in particular VT -symmetric quantum electrodynamics, with imaginary 
electric charge and axial vector potential 17J. 

We already know from previous work on bound states in the Schrodinger 
equation that the VT symmetry of the square well is spontaneously broken [18] , [19] 
and that the perturbative derivation of the C operator below the critical value of 
the imaginary part of the potential is by no means trivial [16] . As far as scatter- 
ing states are concerned, however, we have pointed out in our previous work on 
non relativistic scattering |20| , that, even in the case of spontaneous breakdown 
of VT symmetry, transmission and reflection coefficients for progressive waves, 
travelling from left to right on the real axis, and for regressive waves, travelling 
from right to left, have definite non trivial relations, which are not valid for a 
general complex potential. These relations are connected with the fact that the 
imaginary part of a PT-symmetric potential is an odd function of the space 
coordinate x, so that its integral over the x axis vanishes. 

The present work is thus intended as a phenomenological investigation of 
relativistic features, such as overcriticality, in presence of a spontaneously broken 
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VT symmetry, which would hardly emerge for different reasons, either in the 
case of a general complex potential, or in the case of a potential with exact 
asymptotic VT symmetry |20] . 



2 Formalism 

In the present work we assume a VT -symmetric square well potential 

{0, x < —b (I) 

qiVo-m) , -b<x<0 (II) , 
q(V + iVi) , 0<x<+b(III) {1) 
, x > +b (IV) 

where the real and imaginary depths, Vq and V±, respectively, and the half- 
width, b, are positive numbers , while the elementary charge is assumed to be 
q = — 1 for particles, as the time component of a vector potential in the Dirac 
equation in (1 + 1) dimensions 

i^(x,t)=H D (q)*(x,t). (2) 
Here, the Dirac Hamiltonian, Hd, reads 

d 

H D = V(x)-ia x — +/3m. (3) 

Formula ([3]) is written in natural units, % = c = 1, which will be used 
throughout the present work, and the metric is goo = ~9n = +1- The 2x2 
Hermitian matrices a x and f3 anticommute and are traceless with square unity: 
it is thus possible to identify them with two of the Pauli matrices: the choice 
we make corresponds to the standard Dirac representational 

a x = <J X , = <7 Z , (4) 

particularly suited to the study of the nonrelativistic limit of Eq. @ . 

As is evident from the formulae given above, the solution, ^, to the Dirac 
equation[2]in (1+1) dimensions can be written as a spinor with two components. 
The parity operator, V, and the time reversal operator, T, are to be defined in 
a consistent way. In the adopted Dirac representation, we obtain[22] 

P=e"*P Qt r z , (5) 

where 9-p is an arbitrary constant, and Pq changes x into —x. The Pauli matrix 
a z ensures that the upper and lower components of have opposite parities. 
With formula (O as definition of the parity operator, it is immediate to check 
that fy-p (x, t) = (x, t) is a solution to the Dirac equation ^ with potential 
VV(x)V- 1 = V(-x). 

For the time reversal operator, T, we consistently adopt the following form 

T=e ie ^a z K, (6) 
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where 9q- is a constant and K. performs complex conjugation. W7- (x, t) = 
T\& (x, t) satisfies the equation 

d ( d \ 

- i— T (as, t) = f V* (x) - ia x — +ma z \ T {x, t) . (7) 

For simplicity's sake, we may assume Of — —0-p, so that 

VT =P /C , (8) 

since erf is the identity matrix. Definition ([8]) is consistent with the one com- 
monly adopted in nonrelativistic quantum mechanics (see, e.g., section 4.1 of 
Ref.|2Dl). 

Moreover, it is easy to show that the operator 

C'=e ie c>a y , (9) 
with 9c a real number, meets the condition 

C'H D {q)C'- x = -H D {-q) . (10) 
Therefore, ^>c (x, t) = C'^f (x, t) fulfills a modified Dirac equation 
.Wfc> (x,t) 



Of 



H D (-q)* c >(x,t) . (11) 



Note that C'VT meets condition (|10p . too, provided that the Dirac Hamil- 
tonian, Hd, commutes with VT, i.e. V (x) — V* (— x). 

It is worthwhile to point out that the above definitions are different from 
those commonly adopted in textbooks [10]. In particular, the transformed wave 
function ^c'VT (x, t) = C'VT^ (x, t) satisfies the Dirac equation for "antiparti- 
cles" 

i-q^c-vt {x, t) = H D (-<?) VcvT (x, t) . (12) 

In each of the four regions of the x axis defined by formula (Q]), we search 
for particular solutions, $>(x,t) = $0 {x)e~ lEt , whose spatial part, $0 {%)> can 
be written in the compact form 

,±ifcai _ ( u ± s ) \ n ±ikx 



$o(») = u± (k) ■ e*«* = I ;f £j j • e ±lfc * . (13) 

Direct replacement of formula (fT3"|) in Eqs. (QHSJ) yields for momentum k 

k 2 (x) = (E-V(x)) 2 -m 2 (14) 
and for the ratio, A, of lower and upper components 

^ =± fl Jw+^ 5±A( ' w - (15) 
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where the upper components, u±, turn out to be arbitrary non-zero constants, 
set to 1 for convenience. Adopting the matrix notation of Ref. [21] , the general 
stationary solution, *Bj (x), to the Dirac equation in the J-th region of the x 
axis ( J — I, IV) can be written in the form 

*j(x) = ilj(x)( £ ) , (16) 



j 

where Aj and Bj are constant and 

is the matrix whose columns are the two linearly independent solutions <$>j(x) in 
region J, apart from a normalization factor, which does not affect the derivation 
of transmission and reflection coefficients. Here, Xj = kj/{E — Vj + m) and 
kj = \J{E — Vj) 2 — m 2 . It is worthwhile to notice that A/ = Ajy = A = 
k/ (E + m), with k = y/ E 2 — m 2 > 0, for scattering states, since Vr = Viv = 0, 
while A// = XJjj = A and kn = kjjj = K, since V/j = Vj* n . 

By imposing continuity of the general solution, ty, at the boundary, x^ , 
between regions J and J + 1 



««>-U;:)= n -<w>-U 

it is easy to express the coefficients of the general solution in region 71/ (x — > 
+oo) as linear functions of those in region I (x — > — oo) 

.fe)- MD -(a)- <18) 

where 

m d = nj^ (+b) ■ si m (+b) ■ njlj (o) ■ n H (o) ■ 07/ (-b) ■ (-6) (19) 

is the Dirac matching matrix [2Tj . Since det fij = — 2Aj and det f2y 1 = 1/ det f2j 
are independent of x, it is immediate to check that 

det M D = detfl^-b)/ detQ IV (+b) = Xi/\i V = 1 ■ (20) 

It is worthwhile to mention that, in the matching matrix method applied 
to the 'PT-symmetric square well in the one-dimensional Schrodinger equation 
in a previous work of ours[20] , use was made of a matching matrix, M , which 
expresses the coefficients of the general solution in region I in terms of those 
in region IV; therefore, the M matrix of Ref. [20] is the nonrelativistic limit of 
(AI D ) , easily derivable from formula (119|) . 
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Elementary quantum mechanics immediately yields transmission and reflec- 
tion coefficients in terms of M D matrix elements. For a plane wave travelling 
from left to right (L — > R), we must have Bjy = 0, so that, from formulae 

CM 

_ Aj V _ detM D 1 

Tl - r - ir-^—Mi' (21) 



and 



B I M 2l 

Rl ^ r = M = ~Mg- (22) 



For a plane wave travelling from right to left (R — > L), A/ = 0, so that 



gj _ 1 



and 



Thus, for a PT-symmetric square well, Tl^r = Tr—^, while Rl^r 7^ Rr^l- 
The nonrelativistic case [20: suggests the equality of the two transmission coef- 
ficients as a consequence of the intertwining relation TH = H*T, satisfied by 
the Hamiltonian ([3]) with any local 'PT-symmetric potential. The intertwining 
relation would be broken by a non-local potential[20j 23|; in that case, we would 
have T L ^ R ^ T R ^ L , too. 

Transmission and reflection coefficients are entries of the scattering matrix, 

As a consequence of formulae (|20I2H24|) . S is not unitary. As already ex- 
plained in Ref. 20J, our approach uses PT-symmetric potentials as effective 
potentials, and our purpose is different from that of Ref. [24] . which considers 
them as fundamental and thus searches for a different Hilbert-space metric that 
permits conservation of probability. As shown in the same reference, however, 
one should face, in this latter approach to scattering by localized potentials, 
conceptual problems connected with the non-locality of the metric that ensures 
unitarity of the S matrix. 

From the definition of transmission and reflection coefficients in terms of M D 
matrix elements, it is easy to check that the determinant meets the condition 
|detS| = 1. 

Formulae (|17H19|) allow us to obtain compact expressions of the M D matrix 
elements and, consequently, of the transmission and reflection coefficients (|21t - 
I24[) in terms of real k and A defined in region / and complex K and A defined in 
region //. For the latter quantities it is convenient to use the parameterization 



K = ^^Ze l{v + +v - )/2 , A = v/a+/a_e l(¥, + - v - )/2 , (26) 
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with 

a± = \J(E-qV ±mf + q 2 V?, <p± = arctan (qVi/ (E ~qV ±m)) . (27) 
After some algebra, we obtain 

M# = e- 2tkb | 2 cosh (26 Im if) + (yfy 2 cos (26 Re K) 

+i Im A sinh (26 Im K) ( ^ff ) + Re A sin (26 Re K) ( ) } , 

The M22 matrix element turns out to be the complex conjugate of M-fj 1 

M° = (MB)*, (29) 
and the off-diagonal elements read 



and 



M i2 = i { RC |l|" lA [cos (26 Re K) - cosh (26 Im K)] 
Re A sin (26 Re if) ( ^[jf ) - Im A sinh (26 Im K) ( ^ ) } 



Af£ = » { SsAtoA [ CQS ( 25 Rc jq _ cogh ( 26 j m K y 
Re A sin (26 Re if) ( ^^f ) + Im Asinh (26Imif) (^JXf)} 



(30) 



(31) 



In the Vi — ► limit, corresponding to a real square well, the diagonal matrix 
elements (|28H29p reduce to the corresponding ones of Ref. [2T] , the off-diagonal 
elements pOII31[) differ from those of Ref. [21] by phase factors due to the different 
choice of the origin of the x axis (left edge of the well in Ref. [21], centre of 
the well in the present work): more precisely, M® 2 = e 2lfe& .M^[2T], M 21 — 
e -2ifc6^ / jD gjj^ ag ex pected|20). The square moduli of reflection coefficients are 
obviously not affected by these phase differences. 

It is also of some interest to compute the nonrelativistic limits of the M D ma- 
trix elements, in order to compare them with the corresponding expressions for 
the 'PT-symmetric square well in the Schrodinger equation obtained in Ref. [20] . 
To this aim, we need the limits of the basic quantities kj and Aj for E — ► m + e, 
where e (<< m) is the kinetic energy. In that limit, kj = k 2 IV = k 2 — > 2me, 
X 2 = X 2 V = X 2 - e/(2m), k 2 u = (k 2 IU )* = K 2 -> 2m(e - qV + iqVt), 
X 2 j = (Xjjj)* = (e - qV + iqV 1 )/(2m). Using units 2m = 1, as in Ref. [20], kj 
and Aj coincide. In the same units, one easily verifies, after some algebra, that 
M D NT 1 of Ref. (20], as expected. 

The formulation gives above is suited to the description of scattering states, 
with E < —m or E > +m. The energies of discrete bound states, in the 
— m < E < +m range, appear as poles of the transmission coefficient, Tl^r 
— Tr^l, or, equivalently, as zeros of the M 22 matrix element. For bound states, 
k and A in the asymptotic regions become imaginary, k = ik! = iy/m 2 — E 2 and 
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A = iX' = ik'/(m + E). The equation satisfied by real bound-state energies thus 
reads 

e- 26fe ' (£) { (iw) 2 co S h(2blmK(E)) + (^) 2 cos (2bReK(E)) 
+ Im A(£) sinh (2b Im if (£)) ( ^g^gji' ) 
+ Re A(25) sin (26 Re #(£)) ( ffig^g)!' ) } = 

(32) 

In the 14 — > limit, corresponding to the real square well, formula (|32[) 
reduces to the well-known text-book expression [TO]. 

The scatttering eigenfunctions of the "PT-invariant Hamiltonian (3} are not 
eigenstates of VT: this corresponds to a spontaneous breakdown of VT sym- 
metry. A different scenario would be obtained if the PT-symmetric square well 
potential were the space component of a vector potential: in that case, if one as- 
sumes minimal coupling, p x — > p x + qV x , the time-independent Dirac equation 
reads, with our choice of Dirac matrices 

E^> (x) =H' D ^ (x) = 

By assuming for qV x (x) a square well potential of type ([1]), and repeating 
the calculations of the M matrix as in the previous case, one would obtain the 
following relation between the coefficients of the asymptotic solutions in region 
J (x — » — oo) and IV (x — > +oo) 

( Z ) ='"""" ( B, ) ■ (M) 

For a wave travelling from left to right (L — » R), with the boundary con- 
ditions Ai = 1, B IV = 0, formula (52} yields T L _> R = A IV = e - 2iqVob , 
Rl^r = Bi = 0, while, for a wave travelling from right to left (R — > L), with 
B w = 1, Ai = 0, one gets T R ^ L = Bi = e 2 ^ v " b , R R ^ L = A IV = 0. The square 
well potential thus becomes reflectionless and conserves the probability flux. A 
more general 'PT-symmetric local potential of finite range, V(x) = V* (~x), 
(— i> < x + b < +b), would maintain this property, since, in that case, the argu- 
ment 2qVob of the exponential in formula (|34|) would be replaced with the real 
integral q j\ V R (x)dx, where V R (x) is the real part of V, owing to the fact that 
the imaginary part is an odd function of x and does not contribute to it. 

It is an easy matter to check that the asymptotic wave functions are eigen- 
states of V7~, in keeping with the proof given in Ref. [20] that an exact asymp- 
totic VT symmetry necessarily implies that the potential is reflectionless and 
conserves unitarity. It is worthwhile to note that in this case Tr^l 7^ 
as a consequence of the fact that TH' D ^ -ff^T[2"U], The connection between 
exact asymptotic VT symmetry and potential reflectionlessness was proved 
in Ref.|20j for finite range potentials, but it might be more general, since it 
was already pointed out in Ref. [25] for infinite range potentials of the class 
V{x) =-x 2K + 2 (K = 1,2,3,...). 



dx 



qV x (x) 



<j z m 



, (33) 
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3 Results and Comments 



In order to explore scattering properties of an overcritical 'PT-symmetric square 
well, we have performed several calculations of transmission and reflection co- 
efficients as functions of an increasing imaginary depth, Vi, while keeping real 
depth, Vo > 2m, and half-width, b, fixed. In our calculations, Vq and V\ are ex- 
pressed in unit of particle mass, m, and b in Compton wavelengths Ac = 1/m. 
As an example of our results, Figure 1 gives |Tl_>,r| 2 (= |Tr^l| 2 = \T\ 2 ), 




Figure 1: Transmission and reflection coefficients for PT-symmetric square wells 
with q = —1,6 = 5/m,Vb = 3m. 1st row: V\ = 0; 2nd row: V\ = 0.25m; 3rd 
row: Vi = 0.5m. Al-3 : \T\ 2 ; Bl-3: \R l ^r\ 2 ] Cl-3: \Rr^l\ 2 - 



\Rl^r\ 2 and \Rr^l\ 2 versus energy, E, with m = 1, V a = 3, b = 5 and Vi = 
(top panels), V\ — 0.25 (intermediate panels) and V\ = 0.5 (bottom panels). 
As a general comment, the signature of spontaneous pair production inside the 
well is represented by the transmission resonances at negative energies in the 
— 2 < E/m < — 1 range, adjacent to the bound-state region {—l<E/m< +1), 
where the transmission coefficient may have poles on the real axis, correspond- 
ing to bound states. The half-plane E/m > +1 corresponds to positive-energy 
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scattering, the half-plane E/rn < —4 (= —1 — Vo/m) to negative-energy scat- 
tering. 

The top panels refer to a real well (V\ = 0). In this case, the two reflection 
coefficients are equal : \Rl^r\ 2 — |-Rr^l| 2 = |i?| 2 and unitarity holds: |T| 2 + 
|i?| 2 = 1. The intermediate panels refer to a PT-symmetric well with V\ = 
0.25m : spontaneous pair creation is still sizable, the two reflection coefficients 
differ by more than an order of magnitude at their maxima and unitarity is 
broken. The bottom panels refer to a PT-symmetric well with V\ = 0.5m: 
with increasing V\, spontaneous pair creation is drastically reduced, albeit still 
present; the two reflection coefficients differ by order of magnitudes and \Rl^r\ 2 
is sharply peaked at positive energies (particle reflection), \Rr^l\ 2 at negative 
energies (antiparticle reflection). The potential is neither absorptive (|T^j| 2 + 
\Ri^j\ 2 < 1 at all incident energies), nor generative (|T^j| 2 + \Ri->j\ 2 > 1 at 
all energies), but shows an intermediate behaviour. 

As for bound states, the real well considered in our example has four bound 
states at negative energies and four at positive energies, which still persist when 
the imaginary part is weak, such as V± = 0.25m, considered in the intermediate 
panels of Fig.l. Real bound states begin to disappear at a critical value V\ cr it — 
Q.272m. At V\ = 0.5m bound states with real energies do not exist any more. 
The behaviour of the real bound state spectrum with increasing V± appears to 
be consistent with a spontaneous breakdown of PT symmetry, which deserves, 
in any case, a more detailed formal treatment in relativistic quantum mechanics, 
including a discussion of spectral degeneracies. 

The present work should be considered as a phenomenological exploration of 
properties of PT-symmetric local vector potentials in standard relativistic quan- 
tum mechanics; as such, it might be easily extended to non-local potentials |26j. 
Further developments could involve either the study of the C operator [16j when 
PT symmetry holds , or the search for a more general symmetry than PT for a 
Dirac equation that is not PT-symmetric, as discussed, for instance, in Ref.|27j. 
for the Schrodinger equation. 
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